Edwards measures for powders and glasses 
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Can one construct a thermodynamics for compact, slowly moving powders and grains? A few years 
ago, Edwards proposed a possible step in this direction, raising the fascinating perspective that such 
systems have a statistical mechanics of their own, different from that of Maxwell, Boltzmann and 
Gibbs, allowing us to have some information while still ignoring dynamic details. 
Recent developments in the theory of glasses have come to confirm these ideas within mean-field. In 
order to go beyond, we explicitly generate Edwards' measure in a model that, although schematic, 
is three dimensional. Comparison of the results thus obtained with the irreversible compaction data 
shows very good agreement. The present framework immediately suggests new experimental checks. 



The classical way to go from the microscopic dynam- 
ics to statistical mechanics proceeds in two steps: one 
first identifies a distribution that is left invariant by the 
dynamics (e.g. the microcanonical ensemble), and then 
assumes that this distribution will be reached by the 
system, under suitable conditions of 'ergodicity'. For 
granular systems this approach seems doomed from the 
outset: because energy is lost through internal friction, 
and gained by a non-thermal source such as tapping or 
shearing, the dynamical equations do not leave the micro- 
canonical or any other known ensemble invariant. More- 
over, the compaction dynamics is extremely slow and 
does not approach any stationary state on experimen- 
tal time scales. This raises the question of characterizing 
the typical configurations or the region of phase space 
visited dynamically. 

The proposal of Edwards and collaborators is to 
use an alternative distribution for very gently vibrated or 
sheared granular systems, with the static situation as a 
limiting case. It may be summarized as follows: given a 
certain situation attained dynamically, physical observ- 
ables are obtained by averaging over the usual equilib- 
rium distribution at the corresponding volume, energy, 
etc. but restricting the sum to the 'blocked' configura- 
tions defined as those in which every grain is unable 
to move. This definition leads immediately to an en- 
tropy (in the glass literature a 'complexity') S e d w , given 
by the logarithm of the number of blocked configura- 
tions of given volume, energy, etc., and its corresponding 
density s edw = S e d w /N. Associated with this entropy 
are the state variables such as 'compactivity' X 
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That configurations with low mobility should be rele- 
vant in a jammed situation is rather obvious, the strong 
assumption here is that, apart from the usual statistical 
weights, all blocked configurations are treated as equiva- 
lent — any extra weight of dynamical origin that might 
distinguish them is disregarded. The purpose of this let- 



ter is to argue that this 'flatness ' assumption characteriz- 
ing Edwards' distributions is neither capricious (it leads 
to correct predictions for the compaction dynamics of a 
given class of systems), nor obvious (it does not apply to 
other classes of systems) . To do this we devise a method 
to count the blocked configurations and compute aver- 
ages over them. 

Let us briefly summarize the state of the art. A first 
clue comes from exploiting the analogy between the set- 
tling of grains and powders, as when we gently tap a 
jar with flour to make space for more, and the aging of 
glassy systems j§-||: in both cases, the system remains 
out of equilibrium on all accessible time-scales, and dis- 
plays very slow relaxations. 

In the late eighties, Kirkpatrick et al. |?],|| recognised 
that a class of mean-field models contains, although in 
a rather schematic way, the essentials of glassy phenom- 
ena. When the aging dynamics of these systems was 
solved analytically, a feature that emerged was the exis- 
tence of a temperature Td yn for all the slow modes (cor- 
responding to structural rearrang ements) @Q. For our 
purposes here, T dyn can be defined by comparing the ran- 
dom diffusion and the mobility between two widely sepa- 
rated times t and t w of any particle or tracer in the aging 
glass. Surprisingly, one finds in all cases an Einstein rela- 



tion ((r(t) -r{t w )f) = T dyn 



a(r(t)-r(Q) 



, where r is the 



position of the particle and / is a constant perturbing 
field. While in an equilibrium system the fluctuation- 
dissipation theorem guarantees that the role of Td yn is 
played by the thermodynamic temperature, the appear- 
ance of such a quantity out of equilibrium is by no means 
obvious. Td yn is different from the external temperature, 
but it can be shown to have all other properties defining 
a true temperature . 

As it turned out, despite its very different origin, this 
temperature matches exactly Edwards' ideas: T e d w and 
Td yn happen to coincide for mean-field glass models aging 
in contact with an almost zero temperature bath [0-[l5| . 
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In fact, given the energy E(t) at long times, the value of 
any other macroscopic observable is also given by the flat 
average over all blocked configurations of energy E(t). 
Within the same approximation, one can also treat sys- 
tems that like granulars present a non-linear friction and 
different kinds of energy input, and the conclusions re- 
main the same Jl6| . 

A first partial conclusion is then that Edwards' sce- 
nario is at the very least correct within mean-field 
schemes and for very weak vibration or forcing. The 
problem that remains is to what extent it carries through 
to more realistic models. 

In this direction, there have been recently studies |isj | 
of Lennard-Jones glass formers from the perspective of 
the so-called 'inherent structures' (a partition of the 
phase-space in terms of the blocked configurations |fl7||). 
In this context a 'flat weight' assumption — similar in 
spirit but not quite equivalent to Edwards' — also comes 
into question and is tested in various ways. Though there 
are caveats |T^,|o|, the results are encouraging. 

The path we follow is instead to construct the Ed- 
wards measure explicitly in the case of representative 
(non mean-field) systems, together with the correspond- 
ing entropy and expectation values of observables. We 
thus obtain results that are distinctly different from the 
equilibrium ones, and we can compare both sets with 
those of the irreversible compaction dynamics. 

The first model we consider is the so-called Kob- 
Andersen (KA) model pi]] that, though very schematic, 
reproduces rather well several aspects of glasses 0] and 
of granular compaction p3j ; most important, this model 
is non mean-field. A particle can move to a neighbouring 
empty site, on a three dimensional lattice, only if it has 
less than four neighbours in the initial and in the final 
position. (In these 'hard particle' models the tempera- 
ture is irrelevant, and we set it to one.) The dynamic 
rule guarantees that the equilibrium distribution is triv- 
ially simple since all the configurations of a given density 
are equally probable. However, at densities close to p g 
(~ .88), the particle diffusion becomes extremely slow 
due to the kinetic constraints. In order to mimic a com- 
paction (or aging) process without gravity, we simulate a 
'piston' by freely creating and destroying particles only 
on the topmost layer with a chemical potential p [ p2| . 

i) The dynamic measurements are taken as follows: 
starting from low density, we perform a slow compression 
by raising the chemical potential up to a high value p = 3. 
Since the equilibrium density at p — 3 is much larger than 
the jamming density p g , the system falls out of equilib- 
rium and very slow compaction ensues. We record the 
density p(t) and the spatial structure function gdynifjt) 
defined as the probability that two sites at distance r are 
occupied. We also compute the dynamic temperature 
Tdyn by comparing induced and spontaneous displace- 
ments. This is the set of observables we use for testing 
the different measures, which are obtained independently. 



ii) In the equilibrium measure all configurations 
(whether they are blocked or not) have equal weight. It 
is easy to obtain the exact equilibrium entropy density 
per particle s equ a{p) = —phip — (1 — p) ln(l — p). Since 
T = 1 as mentioned above, ds ">^ ^ = —p. The equi- 
librium structure factor is easily seen to be a constant 
9equii(f) — P 2 '- indeed, one main advantage of this model 
is that it is particularly easy to compare small deviations 
from gequiiif), a notoriously difficult task in glassy sys- 
tems. 

Hi) Finally, we obtain Edwards' measure results as fol- 
lows: we introduce an 'auxiliary model' in which parti- 
cles have energy equal to one if the dynamic rule of the 
original model would allow them to move, and to zero 
otherwise. Performing simulated annealing of the auxil- 
iary model at fixed number of particles is an efficient way 
to sample over the configurations with vanishing frac- 
tion of moving particles the Edwards ensembles struc- 
ture function g e dw{r), and to obtain (Figure |l|) S e d w {p) 
as the logarithm of the number of such configurations by 
thermodynamic integration of the energy of the auxiliary 
model with respect to its temperature. We then compute 

rp-l _ 1 ds cdw {p) 
edw p dp 
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FIG. 1. Gibbs and Edwards entropies per particle of the 
Kob- Andersen model vs. density. At high enough density the 
curves are indistinguishable, and join exactly only at p = 1. 
The slope of the tangent to s e dw{p) for a generic p allows to 
extract T e dm (p) from the relation da 5 dro = ^ — ^-r-r ds '^" 1 . 

vfv dp Ted ro (p) dp 

We are now in a position to compare the long-time 
results of the out of equilibrium dynamics i) with those 
obtained with measures ii) and Hi). Figured shows a plot 
of the mobility X (t, t w ) = ± £a=i £^ =1 *jkffl_^))> 
vs. the mean square displacement B(t,t w ) — 
Jn ELi EfcLi - r fc(^)) 2 ), testing the existence 

of a dynamical temperature Td yn , [ p4[ , in the compaction 
data (N is the number of particles and a runs over 
the spatial dimensions). The agreement between Td yn 
and the Edwards temperature T e d w , obtained from the 
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blocked configurations as in Figure [|, is clearly excellent. 




B (t,tw) 

FIG. 2. Einstein relation in the Kob- Andersen model: plot 
of the mobility x(t,tw) vs. the mean-square displacement 
B(t,t w ). The slope of the full straight line corresponds to 
the equilibrium temperature (T = 1), and the slope of the 
dashed one to Edwards' prescription obtained from figure |l| 
at p(t w ) = 0.848. 
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FIG. 3. Structure functions g(r) — p at density p ~ 0.87 
computed with the equilibrium, Edwards' and dynamical 
measure of the Kob-Andersen model. The three sets of 
data come from independent Monte-Carlo simulations. The 
dynamic structure function (circles) is obtained after slow 
compression raising the chemical potential continuously from 
p = 1 to p = 3 in 10 6 Monte Carlo sweeps. The Edwards 
structure function (open squares) is obtained from the aux- 
iliary model. Although the equilibrium value of g(r) — p 2 
is exactly 0, we also obtain it by a Monte-Carlo simulation 
(full squares) in order to show that the difference in the short 
distance behaviour is not an artifact of the numerical simula- 
tion). The size of the typical error bar on dynamical data is 
shown at r = 3. 

In Figure [|we plot i) the long-time dynamic gdyni^i £)j 
ii) the equilibrium g equ a{r) — p 2 , and Hi) the Edwards 
9edw{i") structure factors. The agreement between i) and 



Hi) is good. 

From the results shown so far, a picture emerges where 
the Edwards measure is able to correctly reproduce the 
sampling of the phase space generated by the out of equi- 
librium dynamics of this non mean-field model. We have 
found, however, that at short times or for excessively 
fast compressions, the quality of the agreement becomes 
worse, possibly due to heterogeneities. We refer to a 
longer, more technical paper for a discussion of these is- 
sues, as well as a study of other models (in particular 
the so-called Tetris model |25j for which one recovers the 
same conclusions as for the KA model) and more techni- 
cal details on our numerical methods |2rf |. 

As already mentioned, Edwards' construction can be 
inappropriate for certain models, even though they may 
have a logarithmically slow dynamics. As a represen- 
tative example of this we consider the low temperature 
domain growth dynamics of a 3D Ising model in a weak 
random magnetic field, a model relevant to many physical 



problems |27|. At large times the domain walls are pinned 
by the field, and the dynamics proceeds by thermal ac- 
tivation. The mean energy decreases slowly towards the 
ground state energy. In a large system, the long-time con- 
figurations are made of domains of 'up' and 'down' spins 
having similar volumes, the global magnetization being 
zero. This is quite different from the equilibrium configu- 
rations at the same energy, which are instead magnetized 
(since the energy is near the ground state energy). 

The question in the present context is therefore 
whether a long-time configuration of (low) energy Eq is 
well reproduced by the typical 'blocked' configuration of 
the same energy. By simulating the corresponding 'aux- 
iliary' model, (with auxiliary energy equal to the num- 
ber of spins not aligned with their local field, i.e. to 
the number of 'mobile' spins), we have checked that this 
is not the case: the blocked configurations consituting 
Edwards' distribution at energy Eq are also magnetised. 
Therefore, neither Gibbs' nor Edwards' distributions de- 
scribe the typical configurations obtained dynamically. 

When is then the flatness assumption characterizing 
Edwards' argument justified? A natural criterion, sug- 
gested by glass theory |2^- |30| , |l^| , consists in studying 
how a system explores its phase space, i.e. its 'chaotic- 
ity' properties. After aging for a time t w , two copies 
(clones) are made of the system, and allowed to evolve 
subsequently with different realisations of the random- 
ness in the updating procedure. We have checked that in 
the KA model the two clones always diverge (the slower 
the larger t w , see figure^), while for the 3D Random Field 
Ising model they do not. It is thus tempting to conjec- 
ture that this form of chaoticity is a necessary condition 
to have flat statistical weights for the blocked configu- 
rations. Note that for this criterion to make sense, it 
should always be applied at non-zero (though weak) tap- 
ping or shearing. The condition of chaoticity is however 
not sufficient: Bouchaud's 'trap model' |9j is chaotic but 
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its fluctuation-dissipation properties are not directly re- 
lated to the density of states. 
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t 

FIG. 4. Mean overlap Qt w (t) between two clones in the 
KA model: the two clones are separated at t w and evolve 
subsequently with different noises. Q tw (t) always decreases 
to zero (the slower the larger t w ), showing that the clones 
always diverge. 

To summarize, our study suggests that the proposal 
made by Edwards does indeed make sense and opens a 
door towards a statistical (thermodynamic) description 
of compact granular matter under very weak driving. In 
order to generalise these ideas to stronger forcing, lower 
chemical potential or higher temperatures (as required to 
analyse the experiments in (3^] ) , one has to learn how to 
go from the concept of 'blocked configuration' to that of 
'metastable state', and this requires other tools The 
inherent structure construction could provide a practical 
shortcut. 

The present setting of the problem immediately sug- 
gests experiments to check these ideas, e.g. by studying 
diffusion and mobility of tracer particles within driven 
granular media. 

Finally, let us note that even in the simplest cases, the 
correspondence between Edwards' distribution and long- 
time dynamics is at best checked but does not follow 
from any principle. The situation is thus as if one would 
have checked that the microcanonical distribution gives 
good results for gases, without knowing Liouville's the- 
orem that proves that such a distribution is indeed left 
invariant by the equations of motion. Such more refined 
arguments would be very welcome. 



ACKNOWLEDGEMENTS: MS is supported by a 
Marie Curie fellowship of the European Commission 
(contract ERBFMBICT983561). This work has also been 
partially supported from the European Network-Fractals 
under contract No. FMRXCT980183. 



[9 



[10 

[11 
[12 



[13 
[14 
[15 
[16 
[17 

[is; 

[19 
[20 



[21 
[22 

[23 

[24 
[25 

[26 

[27 

[28 

[29 
[30 

[31 

[32 



Unite Mixte de Recherche UMR 8627. 



S.F. Edwards, in: Granular Matter: An Interdisciplinary 
Approach, A. Mehta, Ed. (Springer- Verlag, New York, 
1994), and references therein. 

A. Mehta, R.J. Needs, S. Dattagupta, J. Stat. Phys. 68, 
1131 (1992). 

R. Monasson, O. Pouliquen, Physica A 236, 395 (1997). 
See, for example Chapter 7 of: L.C.E. Struik, Physi- 
cal Ageing in Amorphous Polymers and Other materials, 
(Elsevier, Houston, 1978). 

H.M. Jaeger, S.R. Nagel, Science 255, 1523 (1992). 
H.M. Jaeger, S.R. Nagel, and R.P. Behringer, Rev. Mod. 
Phys 68, 1259 (1996). 

T.R. Kirkpatrick, D. Thirumalai, Phys. Rev. B 36, 5388 
(1987). 

T.R. Kirkpatrick, P. Wolynes, Phys. Rev. A 35, 3072 
(1987). 

J.-P. Bouchaud, L.F. Cugliandolo, J. Kurchan and M. 
Mezard, in Spin-glasses and random fields, A. P. Young, 
Ed. (World Scientific, Singapore, 1997). 
L.F. Cugliandolo, J. Kurchan, L. Peliti, Phys. Rev. E 55, 
3898 (1997). 

R. Monasson, Phys. Rev. Lett. 75, 2847 (1995). 

J. Kurchan, in Jamming and Rheology: Constrained Dy- 

namics on Microscopic and Macroscopic Sc ales (1997), 



http://www.itp.ucsb.edu/online/jamming2/, S.F. Ed- 

to be published, 



wards, A. Liu, an d S.R. Nagel Eds. 
cond-mat 9812347. 



Th.M. Nieuwenhuizen, Phys. Rev. E 61, 267 (2000). 
S. Franz, M.A. Virasor o, J. Phys. A 33, 8 91 (2000). 



:ond-mat/9911226 



A. Crisanti, F. Ritort, 
J. Kurchan, J. Phys. Condensed Matter 12, 6611 (2000). 
P.H. Stillinger, T.A. Weber, Science 225, 983 (1984). 
W. Kob, F. Sciortino, P. Tartaglia, Europhys. Lett. 49, 
590 (2000). 

G. Biroli, R. Monasson, Europhys. Lett. 50, 155 (2 000). 
A. Crisanti, F. Ritort, A. Rocco, and M. Sellitto, cond 



mat/0006045 



W. Kob, H.C. Andersen, Phys. Rev. E 48, 4364 (1993). 
J. Kurchan, L. Peliti, M. Sellitto, Europhys. Lett. 39, 365 
(1997). 

M. Sellitto, J.J. Arenzon, cond-mat/0006146, to appear 
in Phys. Rev. E. 

M. Sellitto, Euro. J. Phys. B 4, 135 (1998). 

E. Caglioti, V. Loreto, H.J. Herrmann, and M. Nicodemi, 

Phys. Rev. Lett. 79, 1575 (1997). 

A. Barrat, J. Kurchan, V. Loreto and M. Sellitto, in 
preparation (2000). 

T. Nattermann, in Spin-glasses and random fields, A. P. 
Young, Ed. (World Scientific, Singapore, 1997). 
A. Baldassarri, Tesi di Laurea, Univ. di Roma La 
Sapienza, (1995). 

L.F. Cugliandolo, D.S. Dean, J. Phys. A 28, 4213 (1995). 
A. Barrat, R. Burioni, M. Mezard, J. Phys. A 29, 1311 
(1996). 

E. Nowak, J. Knight, E. Ben-Naim, H. Jaeger and S.R. 

Nagel, Phys. Rev. E 57, 1 971 (1998). 

G. Biroli and J. Kurchan, cond-mat/0005499. 



4 



